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We investigate magnetic excitations in an S = 3/2 Heisenberg model representing two-dimensional 
antiferromagnet Ba2CoGe207. In terahertz absorption experiment of the compound, Goldstone 
mode as well as novel magnetic excitations, conventional magnetic resonance at 2 meV and both 
electric- and magnetic-active excitation at 4 meV, have been observed. By introducing a hard 
uniaxial anisotropy term A(5"')^, three modes can be explained naturally. We also indicate that, via 
the spin-dependent metal-ligand hybridization mechanism, the 4 meV excitation is an electric-active 
mode through the coupling between spin and electric-dipole. Moreover, at 4 meV excitation, an 
interference between magnetic and electric responses emerges as a cross correlated effect. Such cross 
correlation effects explain the non-reciprocal linear directional dichroism observed in Ba2CoGe2 07. 

PACS numbers: 75.80.-|-q, 75.40.Gb, 75.30.Ds, 76.50.+g 



Multiferroic materials have attracted both experimen- 
tal and theoretical interests due to giant magnetoelectric 
effects [iHll- Such strong couplings between magnetism 
and electric polarization (EP) are often realized through 
spin-dependent EPs. For example, in cycloidal magnets 
iT!Mn03 {R = Tb, Dy, and others), EP flops from P|jc to 
P||a by changing a magnetic state from be to ah cycloidal 
state through external magnetic fields [3| • Another exam- 
ple is magnetic resonance induced by oscillating electric 
field, or electromagnon, which is observed in an optical 
spectroscopy at terahertz (THz) frequencies for a vari- 
ety of multiferroics compounds, e.g., PMnOa l5|-l7[ and 
CuFei_xGax02 0. The exchange striction 0,0] and the 
spin current |10l4l2| mechanisms are well known as the 
origins of such spin-dependent EP. 

Spin-dependent metal-ligand hybridization has been 
proposed as an alternative mechanism [l^ [3] ■ EP along 
the bond direction connecting metal and ligand de- 
pends on a spin structure at a metal site S™ in a form 
p,„i oc (Sm • YmiYvmi- At & Spin sitc with no inversion, 
such a mechanism can induce an electric dipole which is 
coupled to the spin, and has a potential to induce novel 
features. In fact, magnetic field dependence of the fer- 
roelectricity observed in Ba2CoGe207 can well be ex- 
plained by introducing this mechanism [isj . 

Ba2CoGc207 is a quasi two-dimensional antiferromag- 
net (Fig. [T] (a)). Below Tjv = 6.7 K, Co magnetic mo- 
ments {S = 3/2) show an antiferromagnetic structure, 
where magnetic moments are aligned in xy-plane due 
to an easy-plane anisotropy In the magnetically 

ordered state, peculiar magnetoelectric behaviors have 
been observed 15|, ll7|. For example, EP along [001] 



shows sinusoidal angular dependence with a period of tt 
for a rotation of the magnetic field B'^'' within xy-plane 
at ^ 1 T. As shown in Ref. [la, such magnetoelec- 



tric behaviors can well be explained by a local electric 
dipole moment which couples to the local spin structure 
of Co atom via the metal-ligand hybridization mechanism 
between Co and O atoms. For a classical spin within 
the xy-plane = {S cosO, S sinO, 0), an EP on a C0O4 
tetrahedron along z is described as = —S^K cos 26, 
which reproduces the experimental results. 

However, there still exist several features to be un- 
derstood. One of them is magnetic excitation property 
observed in an electromagnetic wave (EMW) absorption 
experiment (AE) in the THz frequency regime, which 
indicates magnetic resonances at w ^ jig . 2, and 4 
meV [1^. The lowest two peaks can be assigned to spin 
wave branches which have been reported in the inelas- 
tic neutron scattering experiments (INS) [l3|- Here, two 
distinct modes exist at the F point in the two-sublattice 
ground state due to an anisotropy. However, the ori- 
gin of the excitation at a; ~ 4 meV is not clear within 
magnon pictures. The other point to be understood is 
the THz AE on several EMW polarizations which indi- 
cates that the excitation at 4 mcV is induced by both 
magnetic and electric components of EMW, whereas the 
excitation at 2 meV is excited mainly by the in-plane 
magnetic component. Moreover, the resonance at 4 meV 
shows a non-reciprocal directional dichroism (NDD) un- 
der the external magnetic fields i.e., absorption in- 
tensity strongly depends on the EMW propagation direc- 
tions (forward -f-k or backward — k). In contrast, NDD 
is not clearly observed at the 2 meV resonance. The ori- 
gin of the magnetic excitation and the absorption mech- 
anism is very important to understand the principle of 
the NDD. 

In this Letter, we propose that a uniaxial anisotropy 
term A(S'^)^ (A > 0) gives clues to understand these 
features. In an S* = 3/2 system, the uniaxial anisotropy 
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FIG. 1: (Color online) (a) Crystal structure for Ba2CoGe207. 
Top view of C0O4 tetrahedra is illustrated by squares. An 
5 = 3/2 spin locates on a Co-site with a spin interaction J. 
DM interaction -Df^ (Sf SJ - 5"^ SJ) is also included. On 
each bond, positive (negative) sign of D^j is represented by 
($>>), and the direction from i- to j-site is indicated by 
the arrow. The tetrahedra C0O4 on A- and B-sublattices 
are rotated around z axis with a rotation angle ft and — k, 
respectively. See also Ref. [l^. (b) Magnetization curve for 
B^^'lla;. (c) EP along z-direction in the in-plane magnetic 
field B"" = (B'=''cos<;/)B,-B'"'sin(/)B,0) (B'^" = 5 T) on a 12- 
site cluster. Experimental data in (b) and (c) are extracted 
from Ref. [l^ and Ref. [U, respectively. 



and pI^K [cos(2K,){(S'f )2 - {Sf)^} - sm{2Ki){Sf Sf + 
SfSf )], where is the rotation angle with Ki = k{-'k) on 
the ^(i?)-sublattice as in Fig. [T] (a) Magnetization 
and EP are defined as M'^ = J2 5MsMo5'7 ^.nd = J2Pi 
(7 = x,y, and z), respectively. M and P under the in- 
plane external magnetic field B'^^j] (cos 0b, siiKpB, 0) have 
been calculated by an exact diagonalization on A'^-site 
clusters (A^ = 8, 10, and 1 2) To reproduce A/^ and P" 
observed in Refs. UM and I21I. the parameters in Eq.([T]) 
are estimated as J/A = 0.125, {D'^l/A = 0.005, A = 1.3 
meV, K = 7r/8, and K = 3.8 x 10~^^ C-m. Here we 
use g ~ 2 while V = 1.0 x 10~^® m^ is the volume per 
Co. As typical examples, the magnetization curve along 
B^'^lja; and magnetic field direction dependence of P^ at 
= 5 T on a 12-site cluster are shown in Figs. [T] (b) 
and (c), respectively. Here, system size effects are found 
to be negligibly small. Note that the magnon energy ob- 
served in the INS [3l at 2 K is also reproduced with this 
parameter set (see Fig. [2] (a)). We have confirmed that 
the results are qualitatively robust against choices of the 
parameters within the strong anisotropy limit J/A <C 1. 

Let us consider excitation processes by magnetic com- 
ponents of EMW, i.e.. Ml transitions, which are related 
to the imaginary part of the magnetic susceptibility: 



splits single spin energies into two doubly degenerate 
states: | ± i) with an eigenenergy A/4 and | ± |) with 
9A/4. Here, |m) is a state with 5*2=771 for spin S = 3/2. 
In the strong anisotropy limit J/A <C 1, where we ne- 
glect the higher energy spin states | ± |), an 5 = 3/2 
Heisenberg model can be approximated by an XXZ model 
H<=ft = -£ J{A(crf crj + afa^j) + ala'-} with A = 4 by us- 
ing an 5 = 1/2 spin operator a 120|. It should be a good 
approximation to reproduce the lowest two branches of 
excitations. In fact, A ^ 2.5 gives a good fit to neutron 
data in Ref. [l^ On the other hand, the highest energy 
mode at 4 meV can be assigned to magnetic excitation 
due to the single ion anisotropy gap 2A. 

To clarify the absorption processes at THz frequencies 
in detail, we investigate an 5 = 3/2 Heisenberg model on 
a square lattice with the uniaxial anisotropy term under 
external magnetic field B'^'^: 

H = ^ { J s, . s, - (s^sy - sfs^)} 

n.7i. 

+ ^{A(5n'-.9A'sB-.S,}, (1) 

where is an 5 = 3/2 spin operator on i-site. The direc- 
tions of Dzyaloshinsky-Moriya (DM) interactions Df^ on 
each bond can be determined uniquely from the crys- 
tal structure as in Fig. [1] (a) [l^]. The DM interac- 
tion lifts the two-fold degeneracy in antiferromagnetic 
ordered states. Refiecting the rotation of C0O4 tetrahe- 
dron around z axis, local EPs on i-site are given by p^ = 
-K [cosi2n,)iS^Sf + SfS^) + sin(2f^,)(^f + S^Sf)], 



imx;T(^) = 



hNVfiQ 



^|(n|M^|0)|'^(c.-c.„o).(2) 



Here |0) is the ground state, |n) are excited states and 
hiUno S'le excitation energies to \n) while j = x,y, and z. 
Eq. ^ is calculated on A^-site clusters {N = 8, 10, and 
12) by the Lanczos method [2^, where the (5-function 
is replaced by a Lorentzian with a width e/A = 0.1. 
The results at B'^^ = are shown in Fig. [2] (a). Out of 
plane component Imx™"(<^) vanishes. In-plane compo- 
nents are found to be identical, Imx™(w) = Imx™(w). 
They show that the magnetic components iJ" and Hy 
induce magnetic resonances at around 2 meV and 4 meV 
(Fig. [2] (a)). As shown in the figure, the system size ef- 
fects are small. Hereafter, we show the results on the 
12-site cluster. We indeed see that the excitation around 
2 meV corresponds to one of the spin-wave branches ob- 
served in the INS [l6| , while the higher energy mode is an 
excitation accompanied with the anisotropy gap excita- 
tion 2A. These features are clarified from J dependence 
of the peak positions. As shown in the inset of Fig. [2] (a) , 
in decreasing J, the high energy peak continuously shifts 
to single site gap excitation 2A = 2.6 meV, whereas the 
low energy peak position is proportional to J. 

When spin states couple to electric fields through EP, 
El process may excite magnetic excitations [l^,^]- Such 
processes can be clarified from the dielectric susceptibil- 
ity via spin-dependent EP 
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FIG. 2: (Color online) (a) Im x^Ti^) ^or -y = x, and z. 
Magnon energy observed in INS is extracted from Ref. [TgI . 
Inset; Peak positions around 2 and 4 meV for Im x™"('^) 
a function of J/A on 12-site cluster, (b) Im x'^%{^) for ^ ~ x, 
and z on 12-site cluster. Peak positions observed in THz AE 
for is" and polarizations are extracted from Ref. [Tgl . 



At B^^^ ^ 0, in-plane components of dieleetric suscepti- 
bility are found to be uniform, Imx™(w) = luix^^^^uj), 
as in the case for the magnetic susceptibility. Contri- 
butions to the 2 meV absorption are small. The 4 meV 
resonance is active for any electric components (see Fig. [5] 
(b)). From these results, we conclude that the selection 
rules and the peak positions are consistent with those 
obtained in the THz AE 

The temperature dependence of THz AE can also be 
explained qualitatively. In Ref. [l^ 4 meV absorption is 
observed even above T/v, whereas absorption at 2 meV 
vanishes at T^r upon increasing the temperature. The 
anisotropy gap excitation energy 2A ~ 30 K is larger 
than T^r, and such a resonance can be observed even 
above Neel temperature, i.e., Tjy < T < 2 A. However, 
the resonance at 2 meV vanishes above Tjv , since the spin 
wave excitation exists only in the ordered state. 

In practice. Ml and El processes arc invoked through 
the interaction with EMW as W = -E" • P - H'^ • M, 
where E"^ (H") is the electric (magnetic) component of 
EMW. Provided that both Ml and El processes induce 
an identical excitation, there is a cross correlation be- 
tween magnetic and electric components of EMW, i.e., 
the interference between electric and magnetic responses. 
As we show details in the following, the effects of the in- 
terference can be observed directly as the linear NDD, 
e.g., the interference enhances absorption intensity for 
the EMW with a propagation vector +k but weakens 
that for the EMW with — k, since reversing k is equiva- 
lent to reversing the relative sign of E" and H" due to 
H" = (l/(Uo'^)k X E". As a typical case, wc consider 
dynamical magnctoclcctric susceptibility for and 
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FIG. 3: (Color online) (a) Im xSzi'-^) under the external mag- 
netic fields B^'^lla;. At 7 T, Im x7S{oj) is also shown, (b) 
Difference of the absorption coefficient by EMW propagating 
directions Aa(a;) at B'^ = 7 T, which can be decomposed to 
ferromagnetic component Aap (a;) and antiferromagnetic one 
AaAF(i^). Aa(a;) obtained from THz AE for and H"^ are 
extracted from Ref. [l9l. (c)-(f) Two sublattice spin structures 
in B"^. See text for detail, (g) Excitation process in the limit 
B^P/A « 1. 



|F-|n)(n|A'/-|0)U(c^-c^„o),(4) 



where c = l/-y/eo/^o- The results under B^'^jja; are shown 
in Fig. |3] (a). Iuyx!^^{uj) is enhanced around 4 meV 
excitation. Note that Imx™°(w) is much smaller than 
Im^x^zi^)- IiiiX™''('^) at 7 T is also shown in Fig. [3] (a). 

Experimentally, such a cross correlated effect can 
be observed as the linear NDD [l^. By introduc- 
ing a complex refractive index iV, a polarized plane 
wave with E"||z, H"||a; and k|jy is described as i?" = 
E^eyi^{-iuj{t - (Ny/c))) and i/J exp(-ia;(i - 

(Ny/c))). From the Maxwell's equations, is given 
as v^eirMA^irM ± xS?(w) where fi^xii^) = 

1 + ASr(^) and eM = 1 + Xtli^)- Here, N+ (N-) 
is a complex refractive index for EMW propagating to 
+y (— y) direction. Thus, non-reciprocal part of an ab- 
sorption coefficient a^{uj) = {uj/c)lmN^ is given by 
Aa{uj) = a+ {uj) ~ a- (uj) = (2w/c) ImxS'^ (t^)- Aa(a;) 
at = 7 T is shown in Fig. |3] (b), where the peak po- 
sition and the magnitude of Aa(aj) are consistent with 
those observed in THz AE [Tit . 

Let us consider the magnetic origin for non- 
reciprocal part of the absorption coefficient Aa{Lu) ^ 
{2uj/c)lmx™°{uj). Under the external magnetic field 
B'''^ _L z, spin structure in the Neel ordered state is 
uniquely determined due to an energy gain of the DM 
term, e.g., the state in Fig. |3] (c) is stabilized by B™|la; 
with B^^ > 0. The EMW propagating to — y direction 



4 



in Fig. [3] (c) corresponds to that propagating to +y di- 
rection in Fig. [3] (d) which is realized by a 180° rotation 
around z axis on a spin site. Thus, reversing the magnetic 



field 



is equivalent to reversing the EMW 
direction — )- —A:*', which is consistent with the ex- 
perimental observation [l^. Note that, when we reverse 
the magnetic field — ^ both ferromagnetic mo- 

ment = + and antiferromagnetic component 
= S\ — S^g change their sign as shown in Figs. [3] (c) 
and (d). Each contribution to Aa(w) can be obtained by 
changing the sign of in the calculation, since only m| 
(m^) changes its sign between states in Figs. |3] (c) and 
(e) ((c) and (f)). As a result, Aa(a;) is decomposed into 
two parts Aq;f(w) and AaAF(w), which depend on the 
modification of and m^, respectively. By comparing 
absorptions for the states in Figs [3] (c)-(f), Aq;f(w) and 
AaAF(w) are extracted as in Fig. [3](b). The results in- 
dicate that AaAF(w) is dominant for the NDD around 
4 meV. Generally, NDD can be realized when sponta- 
neous magnetization and EP coexist. In the present 
model, however, realization of a Neel ordered state with- 
out ferromagnetic moment is sufficient to induce NDD. 
Once a single domain structure of the Neel ordered state 
is realized, AQ;(a;) can be finite even at and 

Finally, we note that the selection rules and the cross 
correlated effects can qualitatively be determined within 
a mean field (MF) approximation. The spin Hamilto- 
nian ([T]) can be approximated as TV^^ = '^{^{SiY ~ 
ffMsBf F • S J, where g^BBf ^ = g/isB^ - Y^AH^j) T 
^?w|((5'J>,-(5J),0)} (-(+) for the i-site on A{B)- 
For simplicity, spin states under — 
,5'^'" sine/),, 0) in the limit B^'^/A < 1 are 
Four eigenstates at site i are given in a 



sublatticc) 
{B^^ COS 4 
discussed. 



Iff, 



3 



i)) and |e±) 



form: 

_i^(-g-»3<^,/2|3^_^gi3</,,/2|_3^j<) as in Fig. [3] (g). Eigenener- 

gies for \gf) are A/^T QI^bB^'^ and for both \ef), 9A/4. 
As a typical example, let us consider excitation processes 
induced by if", Hy and E'^ . From the ground state \gi~), 
processes through Sf, Sf and pf are 



Pt\9t) 



COS 4>^\gt) +is\\\(t>,\gl) 
-hA/3(cos0,|e+) -Hisin0,;|e~'))/2, (5) 
s\n(j)i\gf) - i cos (j>i\g~) 
+V3(sin(/),|e+) - icos0,;|e-))/2, (6) 
-V3Kcos{2(j), - 2K,;)|e+) 
-iV3Ksm{2(l)i^2K,)\er). (7) 



We see that i?" can only induce magnetic excitations 
with the anisotropy gap 2A. This is consistent with 
the resuhs in Fig. E] (b). For the EMW with I]'^\\z and 
H" II (cos (j>i^ , sin , 0) , the cross correlation effect is qual- 
ified by a spectral weight A/((/)") cx / (cos (ji^lvax^^ (uj) + 
sin (j>^lmx'^^{uj)) did. The canted Neel ordered state un- 



der B'"'||(coS(^B,sin0B,O) gives 4>, = (jjB T {tt / "2 - M (- 
(-f-) for the spin on A (B) -sublatticc), where 0o is a spin 
canting angle. By applying Eqs. ©-([T]) to Eq. dU, we 
obtain that 



A/ ((/)„) cx cos( 



6b) sin(2K 



(8) 



For 0B = (a state in Fig.[3](c)), we obtain that A/(0) (x 
sin(0o - 2k) (for H"|jx) and A/(7r/2) = (for H"||?/)as 
already expected from the symmetry argument in Ref.ll9l. 
We see A/(0) ^ even for <j>o = 0, which indicates the 
existence of the NDD in a coUinear Neel ordered state. 

In addition, wc can predict the NDD for k||[010], 
E'^IIIOOI], and H"||[100] (0^ = -7r/4) under B°%)10] 
(0B = ''i'/4), although there is no spontaneous EP [l5| for 
this B'^'' direction. In fact, xfioo][ooi] (^) ^d*^!' B'=''||[010] 
is found to be non-zero in the numerical calculation. 
Observation of the NDD in this condition is a crucial 
test for the validity of our theory. As another exam- 
ple, we can easily derive AI^x cx / Imx™''(a;)da; = for 



the EMW with 



and H'" 



under B'^ 



which 



is consistent with the results calculated at B°^ = 7 T: 
/ lmxT.'itu)dtu « / linx7^iLo)di, (see Fig. [3] (a)). 

Our results indicate the potential of the spin- 
dependent metal-ligand hybridization mechanism for 
novel absorption processes which might be observed in 
a wide range of materials with a spin at a site without 
inversion center, e.g., in a tetrahedron and a pyramid of 
ligand atoms. 
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